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Abstract. Aims. We have conducted a study in radio wavelengths and in X-rays of the pulsar wind nebula 
(PWN) in the supernova remnant (SNR) GO. 9+0.1 with the goal of investigating in detail its morphology and to 
accurately determine its characteristic parameters. 

Method. To carry out this research we have observed the PWN at A 3.6 and 6 cm using the Australia Telescope 
Compact Array (ATCA) and combined these data with existing multiconfiguration VLA data and single dish 
observations in order to recover information at all spatial scales. We have also reprocessed VLA archival data at 
A 20 cm. From all these observational data we have produced high-fidelity images at the three radio frequencies 
with angular resolution better than 3". The radio data were compared to X-ray data obtained with Chandra and 
in two different observing runs with XMM-Newton. 

Results. The new observations revealed that the morphology and symmetry suggested by Chandra observations 
(torus and jet-like features) are basically preserved in the radio range in spite of the rich structure observed in 
the radio emission of this PWN, including several arcs, bright knots, extensions and filaments. The reprocessed 
X-ray images show for the first time that the X-ray plasma fills almost the same volume as the radio PWN. 
Notably the X-ray maximum does not coincide with the radio maximum and the neutron star candidate CXOU 
J174722. 8-280915 lies within a small depression in the radio emission. From the new radio data we have refined 
the flux density estimates, obtaining Spwn ~ 1-57 Jy, almost constant between A 3.6 and A 20 cm. For the whole 
SNR (compact core and shell), a flux density S20cm= 11.5 Jy was estimated. Based on the new and the existing 
A 90 cm flux density estimates, we derived a spectral index qpwn = —0.18 ± 0.04 and a? s h c ii = —0.68 ± 0.07. 
From the combination of the radio data with X-ray data, a spectral break is found near v ~ 2.4 x 10 12 Hz. The 
total radio PWN luminosity is L ra( jio = 1.2 x 10 35 erg s _1 when a distance of 8.5 kpc is adopted. By assuming 
equipartition between particle and magnetic energies, we estimate a nebular magnetic field B= 56 fiG. The 
associated particle energy turns out to be U pa rt = 5 x 10 47 erg and the magnetic energy U mag = 2 x 10 47 erg. 
The high ratio between magnetic and particles flux energy density suggests that the pulsar wind just started to 
become particle dominated. Based on an empirical relation between X-ray luminosity and pulsar energy loss rate, 
and the comparison with the calculated total energy, a lower limit of 1100 yr is derived for the age of this PWN. 

Key words. ISMJndividual (GO. 9+0.1) — ISM: supernova remnants — ISM: pulsar wind nebulae — radio contin- 
uum: ISM— X-rays:ISM 



1. Introduction 

Radio composite supernova remnants (SNRs) consist of a 
shell and a spectrally distinct inner nebula, presumably 
a pulsar wind nebula (PWN), powered by the wind of 
relativistic electron/positron pairs from a central pulsar. 
Only in a fe w cases, however, has th e central pulsar been 
detected (see lKaspi fc Helfandll2002l . for a review). 
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Recent observations of several composite SNRs carried 
out with Chandra X-ray Observatory have resolved out 
on arcsec scales complex structures in the interior of sev- 
eral PWNe. These s tructures include torp i ds, axial bipo- 
lar je t s, wisps, etc. (IWeisskopf et al. I l200d: Hclfand et al 
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Roberts et al. 20031 etc.). Images in the differ 



ent spectral domains are essential to understanding the 
physics of PWNe. Particularly, the radio emission depends 
on the history of the nebula and represents the combina- 
tion of the efficiency of the pulsar in providing accelerated 
particles and magnetic fields, and the expansion history. 
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The expansion history, in turn, depends on the density 
and geometry of the medium that confines the relativis- 
tic particles and fields (i.e. the interior of the SNR, that 
includes stellar ejecta and the presence of forward and re- 
verse shocks) . The detailed analysis of the geometry and 
structure of the PWN and the parent SNR, can shed light 
on the coupling mechanisms between the neutron star, the 
relativistic wind nebula and the surrounding SNR plasma. 

GO.9+0.1 (RA= 17 h 47 m 21fi, dec= -28°09', J2000) is 
a composite SNR located in the direction of the Galactic 
center and at about the same distance (assumed through 
this paper to be 8.5 kpc). It is characterized by a bright, 
centrally condensed synchrotron nebula, approximately 2' 
in size, and a weak surrounding radio shell, about 8' in 
size, for which radio spectral indices a coro ~ —0.12 and 
asheii ~ —0.6 (where S oc v a ), have been proposed for 
the core and shell respectively (jHelfand fc Becker! 11987 : 
La Rosa et al~ll2000l ). 

In the X-rays domain, t he first detection was reported 
bv lHelfand fc Beckerl (jl987l ) based on IPC-Einstein obser- 
vations, who concluded that the observed flux could come 
either from the compact core or from a combination of 
core and part of the bright limb of the s hell of GO.9+0.1 



The co re X-ray emission was dete cted bvlMereghetti et aL 
(|l998h using BeppoSAX satellite. ISidoli et all (|2000l ) later 



confirmed these results on the basis of better quality data. 
These early detections are indicative of the presence of a 
young neutron star powering the nebula, although no co- 
he rent pulsations are fou nd . 



Gaensler et al.l (|2001l ) presented the results of 35 ksec 



ACIS Chandra observations of the PWN, between 0.5 and 
8.0 keV. From these images, the authors identify a faint 
semicircular arc and a jet-like feature that define a sym- 
metry axis, which they interpret as evidence of a torus 
of emission in the pulsar's equatorial plane and a jet di- 
rected along the pulsar spin axis. No X-ray emission is de- 
tected in correspondence with the radio-shell nor its inte- 
rior. Based on these observations the authors propose that 
the hard point-like X-ray source CXOU J174722.8-280915 
detected at energies above 3 keV, is the best candidate for 
a central pulsar that would be powering the inner nebula. 
This point source has a rather low ratio of magnetospheric 
pulsar emission to surrounding nebular emission, with a 
luminosity that amounts only 0.5 % of the total PWN 
lu minosity in the energy range 2-10 keV. 

Porquet et all (j2003l ) carried out an X-ray study of the 



PWN within GO.9+0.1 using the XMM-Newton EPIC- 
MOS and EPIC-PN cameras. The images obtained in 
the energy band 1.5-12.0 keV show an amorphous neb- 
ula with a bright maximum towards the cast surrounded 
by extended diffuse emission. At the spatial resolution of 
X MM-Newton (8"), the arc and jet-like features noticed 
bv lGaensler et all (|200l[ ) are not obvious. The X-ray spec- 
trum within the PWN softens from the core to the out- 
skirts, consistent with synchrotron radiation losses of high 
energy electrons as they diffuse through the nebula. The 
XMM-Newton study also reveals sp ectral variation s acros s 
the "arc- like feature" identified bv lGaensler et al. ( 2001 ). 
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Fig.l. Th e SNR GO.9+0.1 at 330 MHz as taken from 
Nord et all (|2004h . 



with the eastern part of the arc having clear indications 
of a very hard photon index (r ~ 1.0), opposite to the 
western part with a very soft spectrum (r ~ 3.2). 

Aharonian et al. I (120051) reported the detection, for the 



first time, of gamma-ray emission in the direction of 
GO.9+0.1 at energies greater than 100 GeV at a level of 
significance of 13 a. The very high energy gamma-rays, 
discovered using the H.E.S.S instrument, appear to orig- 
inate in the pulsar wind nebula. The photon spectrum is 
compatible with a power law with photon index T = 2.4. 

In radio wavele ngths, GO.9+0.1 is promi nent at 57.5 
MHz and 80 MHz dLa Rosa fc Kassiml Tl985T) . It has also 
been observed at 843 MHz (|Gravlll994l ). at 330 MHz as a 
part of the high-resolutio n imaging of the Galactic Center 
regio n dNord et al.ll2004h (Fig [TP. and at 1.5 GHz and 4.8 
GHz (jHelfand fc Beckerll 19871) . 

This paper attempts to analyze the morphology and 
spectral properties of the PWN in GO.9+0.1 in the radio 
range, based on new high-resolution radio images obtained 
from observations at 4.8 GHz (A 6 cm) and at 8.4 GHz (A 
3.6 cm) carried out with the Australia Telescope Compact 
ArrajQ and from reprocessed archival data at 1.4 

GHz (A 20 cm). Also the X-ray emission associated with 
this nebula has been re-analyzed including new XMM- 
Newton data and reprocessed Chandra observations. 



The Australia Telescope Compact Array (/ Parkes tele- 
scope / Mopra telescope / Long Baseline Array) is part of the 
Australia Telescope which is funded by the Commonwealth 
of Australia for operation as a National Facility managed by 
CSIRO. 

2 The VLA of the National Radio Astronomy Observatory 
is a facility of the NSF operated under cooperative agreement 
by the Associated Universities Inc. 
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2. Observations 

2.1. Radio data at 3.6 cm, 6 cm and 20 cm 

The radio continuum emission of GO. 9+0.1 was simultane- 
ously observed at A 6 cm and 3.6 cm using the Australia 
Telescope Compact Array (ATCA) during 12 hours on 
15/16 January 2004. The array was used in the 6B con- 
figuration, which records visibilities from baselines 214 m 
to 6 km. A total bandwidth of 128 MHz split in 32 chan- 
nels was used for each frequency. The absolute flux density 
scale was determined using PKS B 1934-638 as the primary 
amplitude and bandpass calibrator (assuming S6 Cm = 5.83 
Jy and S3. 6cm = 2.84 Jy). Periodic observations of PKS 
B1729-37 were used to correct for changes in gain and 
phase caused by receiver, local oscillator and atmospheric 
instabilities. Image processing at all freq uencies was car- 
ried o ut under miriad software package (jSault fc Killeen 
19991) . 



To improve the uv coverage, the A 6 cm ATCA obser- 
vations were combined in the uv plane with VLA archive 
data acquired in the same radio band with the interfer- 
ometer operating in the hybrid DnC configuration (pro- 
gram AB254 observed July 19 and 22 1984, PI R. Becker). 
Finally, single dish data at A 6 cm acquired with the 
Parkes 64 m telescope (extracted from Parkes-MIT-NRAO 
Southern Survey, Condon et al. 1991) were combined with 
the interferometric data using the MIRIAD task immerge, 
that linearly merges together two images with different 
resolutions after appropriately weighting both data sets 
according to their respective primary beam shapes. This 
addition allowed us to recover all missing flux density due 
to the lack of short spatial frequencies. The faint struc- 
tures associated with the outer SNR shell arc however 
barely detectable because in addition to their intrinsic 
faintness they suffer from attenuation produced near the 
primary beam edge. The astrometry in the final image was 
checked with the 6 cm Catalog of Compac t Radio Sources 
in the Galactic Plane ([White et al.l 120051 ) which has rms 
positional errors lower than 0".67 and 0".84 in RA and 
dec, respectively. The resulting synthesized beam and rms 
noise in the A 6 cm image are listed in Table [1] together 
with the observational parameters for the other analyzed 
wavelengths. It is worthwhile to note that the new A 6 
cm image of GO. 9+0.1 improves in over 50 times the noise 
level with respect to the previously published image at the 
same frequency. 

To produce a high-fidelity image at A 3.6 cm, wc 
searched for more databases at this frequency to im- 
prove the uv coverage. The ATCA data were combined 
with VLA D-array archive data obtained at the same fre- 
quency (observed on March 25 2003, program AJ302, PI 
M. Rupen), after applying an appropriate calibration fac- 
tor. This image is only useful to investigate the central 
PWN since the surrounding radio shell (with a size of ~ 
8') exceeds the primary beam size of the used telescopes 
at this wavelength (~ 5'). 



A new A 20 cm image was produced from the combi- 
nation of archival VLA A-configuration data (observed in 
July 3 1991, program AF209, PI D. Frail) with data ob- 
tained in the VLA hybrid BnC configuration ( observed 
in February 6 and 7 2004, program AY147, PI F. Yusef- 
Zadeh). The A-array observations were carried out at 1465 
and 1515 MHz, using 1331+305 and 1751-253 as primary 
and secondary calibrators, respectively. The CnB obser- 
vations were done at 1385 and 1465 MHz. In this case 
1328+307 and 1748-253 were used as primary flux den- 
sity and secondary phase calibrator, respectively. Because 
of the different observing conditions, both databases were 
separately calibrated and cleaned and later combined us- 
ing the task immerge within miriad software, after decid- 
ing appropriate overlapping ring in the uv plane. In spite 
of the special care put in the cleaning process to mitigate 
the effects of the strong neighbour source Sgr B2 within 
the observed field, some residual striation could not be 
completely removed. 



Table 1. Observational parameters of the radio data 



A 


Instrument 


Beam 


noise 


(cm) 


and configuration 


("x") 


(mJy/b) 


3.6 


ATCA (6B),VLA (D) 


1.5x0.8 


0.05 


6 


ATCA (6B),VLA (DnC),SD 


2.9x1.6 


0.07 


20 


VLA (A, CnB) 


2.5x1.2 


0.12 



2.2. X-ray band 

Wc have analyzed the two existing sets of XMM-Newton 
data on the PWN GO. 9+0.1. The first obse rvation, per- 



forme d in September 2000, was published in lPorquet et al 



(|2003l ). A second longer observation was carried out in 
March 2003. The medium filter was set for all cameras 
(MOS and PN) in the two observations. The data were 
processed using the Science Analysis System (SAS version 
7.1). The periods of high particle background (associated 
with flares) were rejected. Table [5] summarizes the avail- 
able observations and exposure time for each of the EPIC 
cameras, before and after flare screening. Despite the loss 
of a significant fraction of the observing time in the sec- 
ond observation, the statistics is improv ed compared to 
our previous analysis (|Porquet et al.ll2003l ) , allowing a bet- 
ter comparison with the new high spatial resolution radio 
data. 

We also r eprocessed arch i val C handra data, originally 
published by iGaensler et al.l ( 2001 ). to produce an X-ray 
image of the PWN at higher spatial resolution. The obser- 
vation (Obs-id=1036, Seq-num=500102) was performed in 
October 2000 for an exposure time of 35 ks. The image of 
the PWN was produced in the 3-8 keV energy band. The 
tool SMOOTH from the SAS was used to provide an adap- 
tively smoothed image at a signal to noise of 10. 
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Table 2. Summary of the XMM-Newton observations 
used in this paper. 



Obs-Id 


Instrument 


Time 


Time after 
flare screening 


0112970201 


Ml 


17.2 ks 


16.5 ks 


0112970201 


M2 


17.2 ks 


16.2 ks 


0112970201 


PN 


11.8 ks 


9.3 ks 


0144220101 


Ml 


49.4 ks 


25.3 ks 


0144220101 


M2 


49.4 ks 


29.0 ks 


0144220101 


PN 


43.7 ks 


15.7 ks 
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Fig. 2. The SNR GO.9+0.1 at A 20 cm. The angular res- 
olution is 13". 1 x 8". 6 and the rms noise of 0.25 mJy 
beam -1 . 



3. Radio results 

Figure [2] shows the new radio image of the SNR GO.9+0.1 
obtained at A 20 cm based only on the VLA CnB obser- 
vations. The features are similar to those observed at A 
90 cm (Fig. [I]), with the brightest side of the surround- 
ing shell to the west and indication of multiple faint arcs 
in the eastern half. This new sensitive image of GO.9+0.1 
shows considerable diffuse emission in the interior of the 
SNR. 



3.1. Morphology of the PWN 

In Fig. [3] we display the new A 6 cm image of the PWN 
in GO. 9+0.1. Fig. [4] includes the 3.6, 6 and 20 cm images 
of the PWN toge ther with the 90 cm image taken from 
Nord et al.l (|2004h for comparison. 

One of the major contributions of this work is the high- 
fidelity representation attained in the 3.6 and 6 cm images, 
which have revealed that the radio synchrotron emission 
from the confined wind in GO.9+0.1 has a complex mor- 
phology with multiple small and large scale features, in- 
cluding enhancements along filaments, bright knots and 
holes. The various conspicuous features (individualized in 
the image at A 6 cm in Fig. 0| are present at all frequen- 
cies over the wide spectral range between 330 MHz (90 cm) 



and 8.6 GHz (3.6 cm), although the respective brightness 
vary across the spectrum. 



The morphology is dominated by a central band that 
runs approximately from NE to SW mimicking at a larger 
sc ale the "torus" featur e noticed in the Chandra image 
bv iGaensler et al" ( 2001 ). This band, that appears defin- 
ing a symmetry axis, has at least two main maxima (Ml 
and M2 in Fig. [4]), the brightest of which is the western- 
most peak, Ml. From the image at A 6 cm, it looks as 
if this maximum is the result of the overlap of two un- 
resolved concentrations, of which only one (centered near 
17 h 47 m 226.8,-28°09'00") remains visible at A 3 cm. The 
maximum M2, centered at 17 h 47 m 22fi.8,-28°09'00", about 
15" north of the location of CXOU J174722. 8-280915 is 
resolved at A 6 and 3.6 cm, revealing that it is not com- 
pact, but hollow in the center. 



Another interesting feature is the clumpy maximum 
M3, to the south of the nebula. This maximum is ap- 
parently the termination of an almost vertical filament 
that links M2 with M3 (more clearly seen in the images 
at A 90 and 20 cm) and will be discussed below in con- 
nection with the X-ray emission. At the periphery of the 
wind nebula, four striking features can be recognized. The 
northern border terminates in a curious filament pointing 
to the east (named Filament 1 in Fig.U]), while the south- 
ern limb also has a similar narrow filament on the eastern 
corner, but in this case the extension points to the north 
along constant RA ~ 17 h 47 m 256(Filament 2). In addition, 
near the north-western corner a curious small synchrotron 
circular ring is present (Ring 1 in the A 6 cm image, cen- 
tered near 17 h 47 m 18fi, -28°08'40") and close to the south- 
ern border of the nebula but detached from it, a set of 
short filaments can be observed around ~ 17 h 47 m 238 , 
-28°10'30" forming another incomplete ring-like feature 
(Ring 2). The appearance of Filaments 1 and 2, protrud- 
ing from the P WN, resemble the fi laments observed in 
the Crab Nebula ( Hester et a,l.lll996h . and are likely to be 
originated in magnetic Raylcigh- Taylor instabilities at the 
interface between the expanding PWN and its surround- 
ing SNR. 

Our sensitive new images do not reveal any point 
source that could be interpreted as the radio counterpart 
of CXOU J174722.8-280915, the faint hard X-ray source 
proposed as the best candidate for a central pulsar (in- 
dicated by a black cross in Fig. [4j. This is not unex- 
pected however, since a typical young radio pulsar dis- 
tant ~ 8.5 kpc with an assumed luminosity Li.4ghz ~ 56 
mJy kpc 2 (like the median for "high-lumi nosity" young, 
rotat ion-powered pulsars as estimated by ICamilo et al.l 
2006h . would have a flux density Si.4gh z ~ 0.08 mJy, well 
below the sensitivity of the 1.4 GHz image. Moreover, to 
be detectable such a point-like source should be brighter 
than the ~ 70 mJ/beam nebular emission in its vicinity. 
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Fig. 3. Color image of the central PWN in the SNR GO. 9+0.1 at A 6 cm. The angular resolution is 2". 9 x 1".6 and 
the rms noise 0.07 mJy/b. 



3.2. Flux density measurements and spectral study 

In the early pape r repo rting the discovery of GO. 9+0.1, 
Helfand fc Beckerl (|l987l ) estimated the flux density of the 
PWN (called "the core" in their paper) at A 6 and 20 cm 
on the basis of VLA observations. The calculations were 
carried out by integrating the detected emission at both 
wavelengths within a square box 2' on a side and, under 
some assumptions, the authors derive S coro _6cm = 4.16 Jy 
and S coro _20cm = 4.15 Jy. For the shell emission, at 6 cm 
the observational limitations (antennas shadowing limit- 
ing the VLA sensitivity for features on the largest angular 
scales) are compensated by crudely adding a zero-spacing 
flux density to the uv data before mapping. In this way 
the authors derive S s hcii-6cm ~ 7.9 Jy. For the estimates 
at A 20 cm, the authors note that it is not possible to ap- 
ply the same method than before because of the presence 
of the bright source Sgr B2 at the edge of the observed 
field, deriving under some considerations S s hoii-20cm ~ 
14.4 Jy. Finally, these results were combined with data 
collected from the literature between 408 MHz and 22 
GHz to carry out a global spectral study. Here the au- 
thors reasonably note that because of the different beams 



used in the various multifrcquency observations, it is ap- 
parent that some observers measured only the core, i.e. 
the PWN alone, while others included some or all of the 



shell e mission. From this combination, iHelfand fc Becker 
(jl987l ) obtain a corc ~ —0.12 and a s hcii ~ —0.6, concluding 
that GO. 9+0.1 is a composite SNR. 

Later lLa Rosa et al.l (|2000l ) observed GO.9+0.1 at A 90 
cm as a part of the Galactic Center study using the VLA 
(HPBW ~ 43" and rms sensitivity 0.5 mJy beam -1 ). 
In this work the authors report S s hoii-90cm ~ 16 Jy and 
S core -90cm ~ 3.5 Jj[f| concluding, after comparing with 20 
cm VLA data, that a coro ~ +0.12 and a s hcii ~ —0.77, 
which confirms the classification as a composite remnant 
suggested by Helfand & Becker (1987). 

Although it is beyond any doubt that the SNR 
GO. 9+0.1 has two different spectral components, the PWN 
with a flat spectrum and the shell with steeper spectrum, 
it is useful to revise the flux density estimates based on 
the new observations where the contributions at all spatial 



3 Note that in Table 1 the authors list a flux density of 4.8 Jy 
for the core, while in the text the value mentioned for the same 
component is 3.5 Jy, after correcting for diffuse background. 
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Fig. 4. Grey-scale images of the central PWN in the SNR GO. 9+0.1 at 90, 20, 6 and 3.6 cm wavelengths. The black 
cross marks the position of the X-ray point source CXOU J174722. 8-280915. The different features individualized in 
the panel corresponding to the 6 cm image are discussed in the text. 



scales have been adequately considered and processed us- 
ing modern image reconstruction algorithms. Besides, the 
new images with very good spatial resolution allow us to 
accurately determine the areas considered to spatially in- 
tegrate the flux density; also source brightness and back- 
ground contributions are now estimated on the basis of 
very sensitive data. 

In Tabic [3] we list the flux densities estimated from the 
new 3.6, 6 and 20 cm images. The first column lists the 
total SNR emission (that is, including contributions from 
the external arcs, diffuse interior and core), the second col- 
umn corresponds to the core component alone, the third 
one has the same but with the underlying shell emission 
subtracted, and the fourth one lists the shell component 
alone. The errors quoted in Table [3] include the intrinsic 



noise as well as possible uncertainties introduced in the 
selection of the level of background contribution. 

At A 3.6 cm only the flux density of the core compo- 
nent is listed because, as mentioned in §2.1, at this wave- 
length the shell can not be appropriately mapped with a 
single pointing. At A 6 cm, the presence of diffuse emis- 
sion associated with the Galactic plane contaminates the 
SNR flux density estimates. This effect was considered 
by subtracting a background at a level that leaves the 
SNR GO. 9+0.1 emission clearly detached from surround- 
ing emission. Also, only lower limit is listed for the flux 
density at A 6 cm of the whole SNR since, as mentioned 
in Section 2.1, at this wavelength the outer shell contribu- 
tion is probably not completely recovered. Though Tablc[3] 
lists for completeness the core flux density before correct- 



Dubner et al.: Radio observations of the PWN in GO. 9+0.1 



7 




0.8 
0.6 

>. 

3 . 4 

in 

a> 0.2 
_o 

0.0 
-0.2 

1 2 3 4 5 

log v (MHz) 

Fig. 5. Radio spectrum of the PWN in GO.9+0.1. The 
data from this work are s hown by filled circ l es; da ta at 
1400 and 4850 MHz from iHelfand fc Beckerl (|l987t ) and 
data at 330 MHz from La Rosa et alj ( 200C ) are displayed 
as open squares. The best linear fit based on the new data 
plus the 330 MHz (solid line) produces a = -0.22. A 
linear fit to the new data alone (dashed line) produces 
a = —0.14. Data from lHelfand fc Becker (1987) were not 
used in the fitting. An average a = — 0.18±0.04 is adopted. 

ing for diffuse shell emission contribution, all subsequent 
calculations are carried out based on the corrected flux 
density. 

Table 3. Integrated radio flux densities 



A 


SsNR 


SpwN 


ocorr 
D PWN 


Sshcll 


(cm) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


3.6 




1.57±0.40 


1.35±0.50 




6 


>6.90±0.20 


3.20±0.14 


1.45±0.23 


>3.70±0.25 


20 


11.50±0.70 


3.23±0.20 


1.72±0.30 


8.27±0.70 



obtained at three frequencies, we followed different pro- 
cedures to investigate possible spatial spectral variations 
across the PWN. 

First the spectral study was carried out by perform- 
ing the direct comparison of the different images. To as- 
sure that the range of spatial scales measured at each 
frequency was perfectly matched, we applied an appro- 
priate uv tapering and reconstructed the intcrfcromctric 
images. In addition, to avoid positional offsets, the im- 
ages were aligned and interpolated to identical projections 
(field center and pixel size) . We repeated the ratio of the 
images, first at the angular resolution of the considered 
images and second, after degrading the spatial resolution 
to 5"x 5" and to 10" x 10" in order to minimize possi- 
ble bias, like small scale image artifacts, zero level differ- 
ences, background variations, etc., that could mask real 
variation s. Contrary to what w as shown in the X-rays do- 
main by iPorquet et al. ( 2003f ) where the eastern side of 
the PWN was found to have a clearly flatter (and harder) 
spectrum than the southwestern half, we obtained that the 
radio spectral index distribution is practically featureless, 
with no particular morphology/spectrum correspondence 
or tendency. 

To investigate if the lack of conspicuous spectral 
variations originated in the procedure, we repeated the 
study using "tomographic" images (sec for example 



Katz-Stone fc Rudnickl I1997T ). a method where different 



spectral indices are tested and residuals with a spectrum 
flatter (or steeper) than the a tcst are highlighted as darker 
(or lighter) features in an otherwise uniform grey map. 
This method is very sensitive to fine-scale spectral index 
changes. Again the result was that no clear departures 
from a mean spectral index of a ~ —0.2 are evident within 
the PWN. To accurately trace small local variations in the 
spectral index, a set of homogeneous observations acquired 
at different wavelengths using the same instrument and 
with analogous observing conditions, would be necessary. 



Figure [5] shows the global spectrum of the central 
PWN in GO.9+0.1 between 330 MHz and 8400 MHz 
together with the least square fit s to our d ata alone 



2000l 'i 



(dashed line) and to our data plus lLa Rosa et al. 
330 MHz data (solid line). The flux densities published 



bv lHelfand fc Beckerl (|l987l) were not taken into account 
in this new fitting since, for the reasons mentioned above, 
they can be overestimated. In the case of including the 90 
cm data, the PWN spectral index is a PWN = —0.22, while 
based on 20, 6 and 3.6 cm data alone, «pwn = —0.14. In 
what follows we adopt for the PWN apwN = —0.18 + 0.04. 

For the shell component we have compared the 90 
cm image with the new 20 cm image, obtaining a s hoii = 
—0.68 ± 0.07, in agreement with previous estimates. 

The study of spatial variations of the spectrum across 
the PWN is a sensitive tool for understanding the cou- 
pling between the fresh relativistic electrons and magnetic 
fields constantly supplied by the pulsar and the surround- 
ing plasma. Therefore, based on the good quality images 



4. Radio/X-ray comparison 

The improved X-ray images of GO.9+0.1 were used to com- 
pare with the radio brightness distribution. The full com- 
parison with X-ray images, however, lacks either better 
spatial resolution (8" for XMM-Newton) or higher statis- 
tics (only 35 ks of Chandra data) . 

Figure [5] shows the new 3-6 keV XMM-Newton image 
(left panel) and the 3-8 keV Chandra image (right panel) 
compared with the radio image obtained at 6 cm (cen- 
tral panel to facilitate visual comparison with both X-ray 
datasets) . 

It is significant the fact that the maxima in the two 
spectral regimes do not coincide. The X-ray peak is ap- 
proximately 15" southern of the radio maximum M2, 
while at the position of the brightest radio maximum, Ml, 
only weak X-ray emission is detected. 

The improved X-ray images confirm that X-rays 
and radio emitting areas have a comparable extent. 
Good radio/X-ray correspondence can be noticed be- 
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Oim±a3-akeV 




Fig. 6. Reprocessed XMM-Newton and Chandra X-ray images displayed together with the 6 cm radio image. Both 
X-ray images have X-ray contours (white) and a radio contour (black) overlapped to facilitate the identification of 
features and the multispcctral comparison. The location of the pulsar candidate CXOU J174722. 8-280915 is marked 
by a white plus sign. 



tween the term i nation of the "jet-like" feature noticed by 



Gaensler et al.l (j200lh in an approximate north-south di- 
rection and interpreted as a jet directed along the pul- 
sar spin axis, and the radio maximum M3 located exactly 
at the southern extreme of the X-ray feature. Besides, as 
mentioned in §3.1, at A 90 cm and A 20 cm, the radio 
emission strikingly follows the shape and extension of the 
entire X-ray jet. 

Based on the XMM-Newton data we have estimated 
the X-ray flux in : 2.14 x 10 12 erg cm~ 2 s~ 1 between 2 

1 between 4 and 6 kcV, 



and 4 keV, 1.24 x 10 12 erg cm" 2 s 
and 1.56 x 10 12 erg cm~ 2 s~ 1 in the 6 to 10 keV interval. 
Following the procedure used for the radio data, we ap- 
plied a least square fit to the X-ray data, obtaining a spec- 
tral index ax = -1-1, which corresponds to a photon index 
r = 1 — a = 2.1, which equals t he average value betw een 
the hardest spectrum derived bv lPorquet et al. ( 2003 ) on 
the eastern side of the PWN (r = 1.0) and the softest 
one (r = 3.2) on the western part. In Figure [7] we display 
the radio synchrotron spectrum shown in Fig. [5] together 
with the X-ray data. From the intersection of the spectral 
fitting in radio and in the X-ray domain, we can conclude 
that the spectrum must break around = 2.4 x 10 12 Hz. 



5. Energetics of the PWN 

Based on the estimated radio flux density, break frequency 
fb and radio spectral index, we can calculate the radio 
luminosity associated with the PWN in GO. 9+0.1 between 
~ 10 7 Hz and v h . A total L rad i ~ 1-2 x 10 35 erg s _1 , 
is obtained. This radio luminosity can be compared to 



L x ~ 4.7 x 10 a4 erg s" 1 measured bv lPorquet et all (|2003f) 



between 2 and 10 keV, after correcting for d=8.5 kpc. 
The luminosities ratio, a useful parameter because it is 
independent of the distance, turns out to be Lx/ L ra dio ~ 
0.4. 
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Fig. 7. Spectrum of the PWN in GO. 9+0.1 over the range 
10 8 to 10 19 Hz, with a spectral index or = —0.18 in the 
radio regime and ax = —1-1 as derived from least square 
fits to the radio and X-ray data, respectively. The extrap- 
olation of the X-ray spectrum intersects that of the radio 
spectrum at the turnover frequency v\, = 2.4 x 10 12 Hz. 

It is of great interest to determine the energy require- 
ments for the radio synchrotron emission associated with 
the PWN. Energy is stored in the particl es as well as in 
the magnetic field. Following iMoffetl (|l975h we can express 
the total energy of the source as 



U T =U P 



aALB 



-3/2 



VB 2 /8ir 



where L is the total radio luminosity of the nebula, 
V its volume, a represents the ratio between the energy 
in relativistic electrons and in energetic baryons, and is 
assumed to be 1, while A depends on the lower an upper 
cutoff frequencies and can be estimated as follows: 



.4 



d 1/2 (2a + 2) (^ Q+1 /2 _ n,«+V2) 



C 3 (2a + 1) (zV 



a+l 



Dubner et al.: Radio observations of the PWN in GO. 9+0.1 



9 



The total energy as a function of the magnetic field 
B attains its minimum when the magnetic energy is 
approximately equal to the particle energy (equiparti- 
tion situation). It is usually assumed that the magnetic 
field in the source is actually equal to the value which 
gives minimum energy. The minimum total energy is 
U min = 0.50(aAL) 4 / 7 V 3/7 = 7 x 10 47 crg and the mag- 
netic field for which equipartition holds can be calculated 
from: B(U min ) = 2.3(aAL/V) 2/7 = 56 fxG. This value 
ag rees with the mag netic field strength early proposed 
by iHelfand fc Becker! (Il987l ) for this PWN (though they 
assumed a higher t^b), and is comparable t o the value 



estimated for the PWN in the SNR W44 ( Petre et al 



2002ih lAharonian et all (|2005l) estimated a magnetic field 
strength of only 6 fj,G for the PWN in GO.9+0.1, though 
this number is inferred from model fitting to the broad- 
band spectral energy distribution, where the B field is a 
free parameter. 

Based on our calculated magnetic field, the magnetic 
energy is U m ~ 2 x 10 47 erg and the particle energy 
Up ~ 5 x 10 47 erg. The wind magnetization parameter (ra- 
tio between magnetic flux energy density to that in parti- 
cles) is then a ~ 0.4. Such a high value suggests that the 
wind just started to become particle dominated, though 
the transiti on from a > 1 to a < 1 is not yet clear (see 
Arond 1998 . for discussion). 

A rough estimate of the age of the PWN can be es- 
timated from the comparison of the total energy with 
the rate of rotational energy loss of the pulsar. This 
last parameter can be derived from the empirical rela- 
tion between t he X -ray luminosity and B proposed by 
Possenti et all (|2002l ): 



log L x = 1.34 log E - 15.3 

From this relation, an energy loss rate E ~2 x 10 37 erg 
s _1 is derived, thus implying that the pulsar had to be 
continuously injecting energy during, at least, ~ 1100 yrs 
to create the observed PWN. 

As a by-product of the new radio measurements, in 
Figure [8] we show the broadband spectral energy distribu- 
tion (SED) of the PWN in GO.9+0.1 obtained after com- 
bining the radi o data presented in th is paper plus the 330 
MHz data from La Rosa et al. I (l2000h together with the X- 
ray data and VHE HESS measurements (from the public 
HESS database). 



6. Conclusions 

This paper presents new high-resolution and high- 
sensitivity images of the PWN in the SNR GO.9+0.1 ob- 
tained at different radio frequencies. The study is comple- 
mented with reprocessed X-ray images based on XMM- 
Newton and Chandra data. The new radio images have 
revealed interesting structures in the nebula, like bright 
knots, rings and elongated filaments which might be show- 
ing instability regions at the sites where the expanding 
nebula interacts with the surrounding ejecta. From the 
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Fig. 8. Spectral energy distribution of the PWN in 
GO. 9+0.1 from radio to VHE gamma-ray, joining the new 
estimates in the radio and X-ray bands with HESS mea- 
surements. 

comparison of the radio images with the reprocessed X- 
ray images it is found that the X-ray emitting electrons 
largely fills the volume delineated by the radio PWN. Also, 
the new detailed radio images have confirmed the symme- 
try suggested by the Chandra X-ray observations, with a 
bright central band aligned with the X-ray "toroidal" fea- 
ture and a narrow elongated north-south structure that 
appear as the counterpart of the "jet-like" X-ray feature. 
These good radio/X-ray correspondences are, however, ac- 
companied by notable disagreements, the most important 
of which is the separation observed between the radio and 
the X-ray maxima. 

Based on the new radio images, with contributions 
from all spatial scales adequately recovered, we estimated 
the multispectral flux densities and performed a spectral 
study. In Table [4] we summarize these results together 
with other observed and derived characteristic parame- 
ters of the PWN in GO.9+0.1. Our study revealed a quite 
uniform distribution of radio spectral index across the neb- 
ula, with only small fluctuations around the mean value 
of a r = -0.18. 

Table 4. Characteristic parameters of the PWN in 
GO.9+0.1 



S3. 6cm — 
S6cm — 
S20cm = 
Q^radio 
Lradio — 

v h ~ 
B ~ 

Uparticlcs — 
Umagnctic — 

age > 1100 yr 



1.35+0.50 Jy 
1.45+0.23 Jy 
1.72+0.30 Jy 
-0.18 + 0.04 
1.2 x 10 35 erg s" 
2.4 x 10 12 Hz 
56 ^G 
5 x 10 47 erg 
2 x 10 47 erg 



From the combination of observations in the radio 
regime with X-ray data we traced a broadband spectrum 
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which suggests a spectral break at i/j, w 2.4 x 10 12 Hz. 
On the basis of this information, together with the ob- 
served luminosities and the assumption of cquipartition 
between particles energy and Poynting vector energy, we 
investigated the energetics and the magnetic field in the 
nebula. 
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